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We study entanglement and fidelity of a two-qubit system when a noisy holonomic, non-Abelian, 
transformation is applied to one of them. The source of noise we investigate is of two types: one 
due to a stochastic error representing an imprecise control of the fields driving the evolution; the 
other due to an interaction between the two near qubits. The peculiar level structure underlying 
the holonomic operator leads us to introduce the reduced logical entanglement which is the fraction 
of entanglement in the logical space. The comparison between entanglement and fidelity shows how 
they are differently affected by the noise and that, in general, the first is more robust than the latter. 
We find a range of physical parameters for which both fidelity and reduced logical entanglement are 
well preserved. 



PACS numbers: 03.67.Lx 



I. INTRODUCTION 



Entanglement is one of the most striking features of 
quantum mechanics. Its meaning and implications have 
been deeply discussed since the early days of the theory 
but, in the last twenty years, its study has found a new 
impulse due to its application to quantum information 
[l|. Entanglement has a crucial role in quantum tele- 
portation and quantum cryptography and it is believed 
to be essential in quantum algorithms — the system must 
pass through a state of maximum entanglement during 
the computation. Since quantum systems are very sen- 
sitive to perturbations and to the influence of external 
noise, many efforts have been made in order to under- 
stand when this precious resource may be lost and how 
it may be preserved @, H, 0| ■ 

With the aim of a robust and error free manipulation of 
quantum information, geometric quantum computation 
has been proposed about ten years ago [f| @ . According 
to this approach, the quantum information is processed 
by means of operators depending only on global and ge- 
ometric features of the system's evolution. This feature 
has been shown to be powerful in creating and manipulat- 
ing quantum information with high state fidelity [7|. In 
its non-Abelian version, called holonomic quantum com- 
putation, various implementations have been proposed 
for systems such as trapped ions d, p, Josephson junc- 
tion [l(|, Bose-Einstein condensate [llj. neutral atoms 
[l2|], quantum dots controlled by lasers [HI]. The main 
advantage of this approach is robustness of the single 
qubit logical gate against both parametric and environ- 
mental errors. Parametric errors are due to imprecise 
control of the fields driving the system evolution; they 
have been shown to cancel out if the control fields fluc- 
tuate fast enough 0, EE E, Gjl [Hj| ■ Also in the case of 



environmental errors, the holonomic operators show an 
intrinsic robustness H(| • 

In the present paper we extend the holonomic approach 
to a noisy two-qubit system in order to understand if and 
how the entanglement is preserved in holonomic quantum 
computation. The question is not trivial even when, as 
in the case of the present paper, only one of the two 
qubits undergoes a holonomic transformation. In fact, 
it is typical of holonomic systems that the logical space, 
i.e., the space of degenerate states which go through a 
non-Abelian phase change, is embedded in a larger space 
[2l| . In our case the logical space is is embedded in a 
four dimensional space, the system 's Hilbert space, as de- 
scribed in the next section. Since the evolution is not 
confined to the logical space, the entanglement cannot be 
estimated directly in terms of the von-Neumann entropy 
of the subsystems [2gj, rather one should first project 
the evolved quantum state onto the logical space. This 
projection leads to a non-unitary evolution in the logical 
space, that does not, in general, preserve the entangle- 
ment. We also calculate the fidelity and compare it with 
the entanglement, when the system is affected by noise. 

We consider two types of noise: a parametric error that 
goes off whenever the driving fields are turned off, and 
a "coupling error" due to undesired interactions between 
the two qubits. We show that the holonomic operators 
are robust under parametric error. Moreover, we show 
that entanglement and fidelity are differently affected by 
these noises and that entanglement is more robust than 
fidelity. As one may expect, entanglement is influenced 
by the coupling between the qubits, however, we are able 
to single out a regime in which both entanglement and 
fidelity are preserved. 

The paper is organized as follows. In Sec. [IT] we give 
a brief review of the holonomic approach to quantum 
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computation; in Sec. IIIII wc define the reduced logical 
entanglement and we calculate it in the presence of a 
parametric noise. In Sec. IIVI we evaluate the fidelity. 
In Sec. |V| we analyze both entanglement and fidelity in 
the presence of an undesired coupling between the two 
qubits. In Sec. [VI] we conclude. 



II. HOLONOMIC COMPUTATION 

We briefly review the main features of holonomic quan- 
tum computation for a system (hereafter called system 
A) described by a four-dimensional Hilbert space Ha, 
with a level structure of three excited degenerate states 
\i) (i = 0, 1, a) at energy e, and ground state |G), set for 
convenience at energy 0. 

The system is driven by time-dependent laser fields, 
with AC frequency in resonance with e, inducing transi- 
tions between ground and excited states. In the rotating 
frame representation the Hamiltonian governing the evo- 
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lution of the system is (ft = 1) [9[ 

H (t)= mt)\i)(G\ + h.c.) , (1) 



where f2i(i) are the time-dependent Rabi frequencies of 
the laser fields. The Hamiltonian Ho(t) has four eigen- 
states: two bright states 



\Bi{t)) = j-^mm+EiSum) 



(2) 



where 



n(t) 



and two dark states 



t=0,l,a 



(3) 



\Di(t)) 
\D 2 (t)) 



n„(t)(ni(t)|i) + n (t)\o)) - (n(t) 2 - \n a (t)\ 2 )\a) 



n (t)|i) -^i(t)|o) 
v/|ni(t)| 2 + |fioWI 2 



\n (tw 



(4) 



The bright states have energy E±(t) = ±f2(t) (the posi- 
tive value is associated with |-Bi)), while the dark states 
have zero energy. The time evolution operator associated 
with |T|) is 



U (t) =Te- i ^ Tff °W 



(5) 



where T is the time-ordering operator. 

The modulation of the phase and the intensity of the 
laser fields drives adiabatically the Hamiltonian along a 
loop in the space of the Rabi frequencies — the parameter 
space — with Hq(0) = Hq(T), where T is the final adia- 
batic time. We shall assume that the Rabi frequencies 
have the following time-dependence 



n (t) = 0. sin 6 (t) cos <f>(t) 
= fl sin 9 (t) sin (j)(t) 
Q a (t) = flcos6{t) 



(6) 



where is a constant. This means that we assume that 
the parameter space is the surface of a sphere of radius f2 
and that the adiabatic loop represented by ([6]) is a closed 
curve on it. Note that the energies E±(t) of the bright 
states do dot depend on time and have constant value Q. 
We shall assume that t = and t = T correspond to the 



north pole on the sphere. Thus, from @, and ©, 

|ZM0)> = |£>i(T)) = |0), (7) 
\D 2 (0)) = \D 2 (T)) = \1), (8) 

|Bi(0)> = \B 1 (T)) = -L(\a) + \G)), (9) 

\B 2 (0)) = \B 2 (T)) = -L(\a)-\G)). (10) 

Under the adiabatic condition OT 3> 1, the adiabatic 
theorem guarantees that any superposition of the dark 
states at t = will end up in another superposition of 
dark states at t = T, and that this transformation is 
realized by the unitary holonomic operator 0] 



(11) 



obtained from ([5]) in the adiabatic limit OT 3> 1; A(t) in 
([lip is the connection operator with matrix elements 
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l3 (t) = (D i (t)\j t \D j (t)) 



i,i = 1,2. 



(12) 



Equation (fTTj) is the core of the holonomic approach to 
quantum computation: in view of (|7 |) . (|8 |) and pip , the 
states |0) and |1) can be regarded as logical states span- 
ning the two-dimensional logical space Ca C TLa (the 
"A-qubit") on which Uq(T) acts as logical operator. 
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FIG. 1: (color online) Left: Loop on the sphere of radius f2 in the parameter space which produces the desired logical operator 
with rj = 7r/2. In the calculation of the integral (|15|l the loop begins in (6 m , 0) (with 9 m near to the north pole), the evolution 
is along a meridian up to (6m, 0), then along a parallel up to (#m,0m), along a meridian again up to (0m,<?W) and finally 
along a parallel back to (0 m , 0). We recover the desired loop, that starts and ends at the north pole, in the limit 9 m — > 0. The 
angle variables are supposed to depend linearly on time. Right: Entanglement as a function of ojq, for a = 0, with wi = 0.4 
(dot-dashed), 0.6 (dashed), 1 (solid line). 



The main feature of the holonomic approach consists 
in the fact that Uq(T) depends only on the area spanned 
by the curve ©. To see how this comes about, let us 
evaluate the RHS of (fTTj) : for the matrix (fT2)) one finds 



A(t) 





<p(t) cos 9(t) 



-<j>{t) cos 0(t) 




whence, 



where 



Uo(T) 



n 



cos rj sin r\ 
— sin i] cos 7/ 



dt(j)(t) cos0(t) 



(13) 



(14) 



(15) 



is the solid angle spanned on the sphere during the evolu- 
tion. Therefore, if one wishes to construct, for example, a 
logical NOT operation, one has just to drive the external 
fields in such a way that the closed curve ([6]) spans a solid 
angle rj = it/2. In fact, for this value of rj , the matrix 
(fT4"|) does what a NOT should do, namely, it exchanges 
the logical qubits |0) and |1) (modulo a sign). 

In order to evaluate the integral (p~5|) . we proceed in 
the following way: 

1. We choose a curve as in Fig. []] consisting in a se- 
quence of evolutions along meridians and parallels 
where 9m, <Pm arc the maximal angles spanned 
during the evolution. Therefore, the solid angle 
spanned by the curve is rj = 4>m(1 — cos 9m)- 



2. Since the north pole of the sphere is a singular point 
of the parametrization (JSJ), we choose a loop along 
meridians and parallels which starts (and ends) at 
the angles (£? m ,0), (Fig. [J). The value of the inte- 
gral for the loop starting and ending at the north 
pole, is recovered in the limit 9 m — > (23[. 



We shall now include the effect of parametric noise 
in the model. In view of the adiabatic evolution, it is 
reasonable to assume that the lasers are stable up to 
a relative error in the ratio between noise and signal. 
Thus the effect of an imprecise control of the driving 
fields is easily modeled by replacing the Rabi frequencies 
fii(i) in © with fii(i) = Qi(t) + 5tii(t), where 6tti(t) 
are random fluctuations proportional to the intensities 
fij(rj) with SCli(t) <C Oj(t). (From now on, for easiness of 
notation, we shall omit the explicit time dependence in 
the formulas and write, e.g., 6Qi instead of <5Qj(t).) 

To compute the perturbed geometric operator U(T) by 
means of the perturbed connection A, we pass from the 
Cartesian coordinates J7j + to the spherical ones, 



fi = n + An, 

9 = 9 + A9, 

4> = (j} + A(j). (16) 
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From ©, 



n = l^cii + sfti) 2 , 



arccos 



arctan 



(17) 



By straightforward computation, one obtains the series 
expansion for AQ, AO and A(f> in the small parameters 



An = sn + s 2 n + o(s 3 ) , 

A6 = 69 + 5 2 6 + 0(S 3 ), 
A<j) = 5<j) + 5 2 (l) + 0{5 3 ) . 



(18) 



Here the "<5 n " -terms collect terms of order n in Sili, e.g. 

sn n n a + 50i0 1 n a + sn a (n 2 a - n 2 ) 
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The perturbed solid angle is written as 

?7=?7 + <5?7 + <S 2 ?/ + 0(<5 3 ), 
with first order correction 



(19) 



S r\ 



dt(<f)89 sin + 5<j) cos ( 



(20) 



The perturbed dark states \Di) are easily obtained 
from (dJ) by replacing flj with fl t . Up to the second order 
in 5Qi and at the final time t = T, one finds 

\Di) = (l-^-^)|^) 

+ {5<j> + S 2 ^) - (IB,) + \B 2 )) , 



V2 



\Do 



• + 6 2 cj ) )\D 1 ) + (l 



5 2 A. 



(21) 



Note that \Di) is a superposition of both unperturbed 
dark and bright states thus leading to a population leak- 
age from the unperturbed dark space to the unperturbed 
bright space. 

The perturbed final operator U (T) , written in the basis 
{|-Di), I-D2)}, can be read directly from equation (fT4|) . by 
replacing 77 with rj = r\ + br\. In particular, for 77 = 7r/2, 
up to second order, 



U(T) = 



-5rj — S 2 i] 1 



(Sv) 2 



(22) 



-5f] — 5 2 rj 



III. ENTANGLEMENT IN MULTILEVEL 
SYSTEMS 

We shall now analyze a two-qubit system. We begin 
by specifying the model that we shall investigate: 

1. We consider the composite system formed by two 
replicas of the system considered in the previous 
section — one will be called system A and the other 
system B. The Hilbert space of the composite is 
then Ha <8> Hb, where Hb, the Hilbert space of 
system B, is a replica of TLa- 

2. The logical space of the composite is La ® C-b, 
where, as before, La is the logical space of system 
A, spanned by the logical states |0) and |1), and Lb 
is the replica of La contained in TIb ■ Following the 
standard terminology, we shall call La the A-qubit 
and Lb the B-qubit. 

3. We shall assume that only the A-qubit undergoes 
information processing (for sake of concreteness, we 
shall assume that the A-qubit performs a logical 
NOT operation). 

4. The initial state of the system is taken to be the 
maximally entangled state 



* = i=(|00> + |ll» 



(23) 



We note that the relevant simplifying assumption is 3, 
while 4 is rather standard for the kind of problem we 
wish to address here (see below). Note that assumption 
3 allows two possibilities for representing the time evolu- 
tion of the two-qubit system: 

a) a unitary transformation of the form U <S> I, where 
U governs the dynamics of the system A (e.g., as 
given by (|2"2"jl ). and / is the identity in Hb', 



b) a unitary transformation which accounts for the in- 
teraction between the two qubits. 

Case a) will be addressed in this section and case b) in 
Sect. El 

Our aim is to study how the entanglement of the two 
qubits is preserved by the dynamics and to do this we 
need a quantitative estimate of the entanglement. 

When a system is composed by two sub-systems A and 
B, the entanglement of the composite system, in the case 
of pure states, can be estimated in terms of the von Neu- 
mann entropy 



£ = -Tr (p A log 2 p A ) = -Tr(p B log 2 p B ) 



(24) 



where pa an d pp are the reduced density matrices of the 
sub-systems [22j . However, since in our model the qubits 
are embedded in a larger space, formula (|24p cannot be 
directly applied. 
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To arrive at a suitable estimator of entanglement, we 
proceed as follows: Let P denote the projector on the 
logical subspacc La of Ha and let P®P be the projector 
onto Ca® Cb- Consider 



and 

Then the quantity 



p <g> PpP <E> p 



PA = Tv Bp r ■ 



-Tr(p A log 2 p A ). 



(25) 



(26) 



(27) 



is analogous to (|24|) , and can be regarded as an estimator 
of the fraction of entanglement in the logical subspace 
(see also [24[). Hereafter we shall refer to to £ r as the 
reduced logical entanglement or, when no ambiguity will 
arise, simply as the entanglement. 

We wish now to obtain a convenient formula for £ r 
(under the assumption a) specified above). Firstly, we 
evaluate 

p\ = Tt b {P ® PU ® I\^){^\U* ® IP ® P) (28) 

without relying on any specific form of the unitary U. To 
this end we note that a generic transformation U of the 
logical qubit basis can be written as 



U\0) = PU\Q) +P ± U\0) 
U\l) = PU\1) +P ± U\1) 



(29) 



Let 



/u = \\PU\0)\\ , y/Zi = \\PU\1)\\ (30) 

and denote by \tpo), \ipi), IV'o"}; and IV^i') the vectors 
obtained by normalization from PU\Q), PU\l), P ± U\0), 
and P ± U\1) respectively. Then 



u\o) = v^lv»o) + vi^m > 

u\i) = + VT=Z*W>£) ■ (3i) 

It is useful to introduce the scalar product 

a = (V^ol^i) ■ (32) 



Note that the coefficients VI — represent the leakage 
of the populations and that in general a is different from 
zero. Therefore 

1 



and 
PU®I\V) 



+V^IV>ii} + \/T^i>ri>] , (33) 



y/ujo + oji 



Noting that \ipi) can be decomposed into its components 
along and orthogonal to |-0o) : IV'i) = ctlV'o) + 0\tfo) with 
(ipo\ipo) = 0, a given in ([32]) . and (3 = \/l — a 2 , we find 

/t^l^oO) + a^/wTl^ol) + (iy/uJl\i! l) 



pu®i\v) 



y/u + LOl 



(35) 



Thus, we may evaluate ([2"5]) in the {|^o) > IV'o)} basis, 



UJlCt + OJq 



Pa = 



whence 



where 



W ° + Wl V a^l - - (a 2 - 1) wi 



f' = -(A_log 2 A_+A + log 2 A+), 



(36) 



(37) 



w + wi ± v/w 2 + 2 (2a 2 - 1) ujicjq + lo\ 

A ± = 7Tf — ; — \ (38) 

2 (uiq + uji ) 

are the eigenvalues of p r A . 

In Fig. Q] (right), £ r is plotted as a function of ujq for 
a = and different values of u>i. It follows immediately 
from (|37|) and (f38|) that £ r = I for a = and cj = 
(this includes the trivial case of no leakage for which ujq = 
Wi = 1, as well as non-trivial cases of population leakage 
with ojq = uj\ < 1). For a ^ 0, it is always £ r < 1. 

We now specialize to the case J7 = U(T), with E/(T) 
given by (f22|) . Keeping the contributions up to the second 
order in the perturbations SQi, we obtain 

PU(T)\0) = ^T \i> ) = {5^ + 5 2 4>-5r l -5 2 r,)\{)) 



P^(T)|0) = VT^|^> = ^ (|o) + |G)) 



PU(T)\1) 



wi\ipi 



5<j) 2 ((5r/) 2 



•<ty<$7?)|0) 



+ (64> + 5 2 4>-5r 1 -5 2 ?7)|1) 



p^c/(T)|i) = yr - ^!^) 



<56» + 5 2 6» 



and 



V2 

a = 0(S 3 ) 
uj = l + 0(<5 4 ) 



(\a) + \G)) 



(39) 



0(£ 3 ). 



(40) 



In the above equations 80 and 6(f> are evaluated at < = T. 
From and (fT9|) it follows that 50 = and that one 
can set 8<p = 0. Thus a = and u>o = u)\ = 1, whence, 
as above, £ r = 1. 
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IV. FIDELITY 

If the system is subject to noise sources, the fidelity 
T is used to quantify the "distance" between the per- 
turbed final state W(T)^, due to a perturbed evolution 
U, and the final unperturbed state Uq(T) ® I*S> , with 
Uq{T) given by tH}. For composite systems, T provides 
an estimation of the performance of logical operators that 
supplement the information already provided by £ r . The 
explicit formula for the fidelity is 



F=\(U(T)i&\Uo(T)®m}\ 



(41) 



For the case we wish to consider first, U (T) = U(T) <X> /, 
with U (T) given by (|2"2")h From (|3"9")l . by a straightforward 
computation, one obtains 



{sef {5<t>) 2 [8nf 



+ 54>5i] + 0{5 3 ) . (42) 



Since, the terms 86 and 8<p are evaluated at time T, their 
contribution is zero and thus 



T ~ 1 



(43) 



The fidelity depends on the error 8r] relative to the area 
spanned during the perturbed evolution, it then appears 
to be less robust than entanglement. 

We shall now estimate 8r] for a simple model of noise. 
We first integrate by parts the second term in Eq. (j2"0|) 
neglecting the contribution at the endpoints 



8rj 



dt(<j>86 + 65(f)) sin 6. 



(44) 



We consider the loop in Fig.[T]passing near the north pole 
and then take the limit for 6 m — ► 0, in the notation of 
the figure caption. We separate the contribution of four 
different parts depending on the evolution along meridi- 
ans and parallels (cf. Fig. [TJ and denote by J i f(t)dt the 
contribution along the ith path. We have 



8r\ 



sin t>M 



1,3 



sin ( 



dt<t>86 . 



(45) 



Note that in the limit of 8 m — » the fourth integral gives 
no contribution. We now suppose that the perturbations 
6Qi fluctuate randomly over a time scale r <C T and go 
off whenever the driving fields are turned off, i.e. that 
they are of the form 



k=0 



k > 



(46) 



where h(x) is the "box" function, which is equal to 1 
in the interval < x < 1 and zero elsewhere, n = T/t, 



and Zl are independent Gaussian random variables, with 
zero average and variance a 2 . 

Using ^ and (HHJ) we can explicitly write the contribu- 
tions for the separate evolutions. Along the first meridian 
we have no contribution since = 0. The only non-zero 
contributions are along the second and third paths; we 
suppose that 9 and <p depend linearly on time (see also 
[23[) and write 9 = ve and (f> = —v^. We thus obtain 

5rj = Vef, cos 9m sin 2 6m J dt (Zq cos 2 + Z\ sin 2 — Z a ) -+ 

- ^ sin2c/> M / dt {Z Q - Z x ) sin0 . (47) 

Inserting (|4"6")l in (|4"T|) one can perform the integration 
and the expansion at first order in t/T = 1/n <C 1 

Srj = ictjM cos 9m sin 2 9m x 



n-i 

k=0 

2 



Z% cos 2 {krv^) Zl sin 2 {krv^) Z% 



M 



)sin(20 M )5Z 



sm(krve)(Z° k -Zl) 



k=0 



Using the independence of the random variables Zl , 
the central limit theorem it follows that for n> 1 



2 ? 
IK ^2 _ 0~ T 

(Sri) w — = — , 
n 1 



(48) 



from 



(49) 



which expresses the cancellation effect already discussed 
in [TH for large n. Thus the fidelity (|4"3"| becomes 



T k, 1 



2 T 



(50) 



The above calculation is valid for any loop in the pa- 
rameter space which moves along a parallel and meridian 
and spans a 7r/2 solid angle. Among these loops, there 
is one which is particularly interesting. If 9m = tt/2 and 
<t>M = tt/2, we are moving on the fl a — &o plane, along 
the equator and on the £l a — Oi plane. As can be seen 
from (|4"T|) , with this loop the contributions along the sec- 
ond and third paths are zero and Srj = independently of 
the characteristics of the noise (variance and correlation 
time). In other words, this loop makes the system com- 
pletely robust against this particular perturbation. This 
fact has a simple geometrical and physical interpretation. 
We recall that the loop describes the way we turn on and 
off the lasers. The above loop is the one in which we 
have always one laser completely turned off while mod- 
ulating the intensities of the other two; for example, in 
Fig. [2] is shown the evolution along the equator (path 1) 
when f2 a laser turned off and we are modulating the f^o 
and intensities. Along this path, Sfl a — while we 
have perturbation of the other two lasers. However, the 
perturbations SQq and Sfl± produce only radial and 89 
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path 1 



FIG. 2: Evolution along the equator (path 1 with 6m = t/2) 
and along a generic parallel (path 2). In the first case, only 
the perturbation Silo and <5S!i for the laser turned on are 
present. 



perturbation which do not affect the solid angle along 
this path. In other words, choosing this particular loop, 
along the single paths, we eliminate the part of pertur- 
bation which can modify the solid angle spanned and as 
result the angle is spanned without error. 



V. COUPLED QUBITS 

In a more realistic situation, the two qubits can also 
interact. This interaction allows to manipulate the two 
systems as a whole and it is the basis for constructing 
two qubit gates. Ideally, one can control the coupling 
strength and turn it on and off depending on the logi- 
cal gate. However, if the interaction cannot be perfectly 
turned off, its presence results in a new source of noise. 
Here, we choose to describe it with a simple model with 



ffj=x|ll><ll| 



(51) 



As before, only the A qubit undergoes an information 
process. Due to the specific form of the coupling (|5"Tj) 
only the state |1) of qubit B will feel the additional in- 
teraction. It is then convenient to analyze the evolution 
in the dark-bright basis for the qubit A and in the state 
|1) for the quit B: {|Z>il}, \D 2 1), \Bil), \B 2 1)} (neglect- 
ing, for the moment, additional errors induced by the 
imprecise control). If fi x, an d the system starts in a 
superposition of dark states, the transition to the bright 
states are negligible and the evolution stay in the dark 
space {|Z?il), | Z?2 1) }• The Hamiltonian Hj in this basis 
is 



H, 



cos 2 8 sin 2 - 



cos W cos sin c 



cos a cos sm ( 



cos 



(52) 



Thus, one sees that Hi breaks the energy degeneracy 
Ei = E 2 = of the two dark states \Di) and \D 2 ), 
relative to the unperturbed Hamiltonian (fT]). The new 
eigenvalues are easily evaluated by diagonalization in the 
presence of Hi , 



Ei = , E 2 = x(cos 2 + 



cos 2 6 sin 2 i 



(53) 



The corresponding eigenvectors can be written as a su- 
perposition of the unperturbed dark states 



with 



\Di{t)) = a(t)\Di(t)} - (3(t)\D2(t)) 
\D 2 (t)) = p(t)\Di(t))+a(t)\D 2 (t)}, 

cos 4>(t) 

y^cos 2 <j){t) + cos 2 0(t) sin 2 <j)(t) 
cos 0(t) sin <j)(t) 

y^cos 2 <j){t) + cos 2 9(t) sin 2 <j)(t) 



(54) 



a(t) 
Pit) 



(55) 



These coefficients satisfy boundaries conditions, related 
to the closed loop in Fig. © with 0(0) = 6{T) = 0(0) = 
<KT) = 

a(0) = a(T) = 1 , /3(0) = @(T) = , (56) 
and normalization 



a 2 {t)+(3 2 (t) = 1. 



(57) 



Note that the energy shift (|53|) . induced by the pertur- 
bation, may produce important modifications during the 
time evolution since non-Abclian effects are based on the 
assumption of degeneracy. If %T 3> 1 the perturbed dark 
states are separated during the evolution and no holo- 
nomic operator (jlip can be produced. For this reason we 
shall focus, in the following, on the more relevant case 
\T < 1, which preserve the logical space. 
The evolution operator is now 



U X (T) = Te~ i ^o[ Hl ^- iA ^] dt 



(58) 



(Note that the first term in the exponential, which rep- 
resents the dynamic contribution of the perturbed dark 
states, is absent in (fTTj) since there the dark states are 
degenerate). Using the new basis {|l)(i) 1 l), \D(t) 2 l)}, 



A ij (t) = (D i (t)l\j t \D j (t)l) i,j = 1,2, 



and thus 



A(t) 



act 

a/3- 



■f3(3 
Pa- 



COS! 



COS f 



(3a — a(3 ■ 
aa + /3$ 

(59) 

The integral L A(t)dt in (|5^|) can be evaluated using 
([S"5)) , (|56|) and (|57|) . The integrals of the diagonal terms in 
(|59[) are always zero because they represent the derivative 
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of the norm of the dark states, which is time independent. 
For the off-diagonal part, it can be easily shown that 
a/3 — (3a — a/ (I — a 2 ) 1 / 2 gives no contribution once it is 
integrated along the closed loop. Thus 



We shall consider as initial state the maximally entan- 
gled state (Hg. From Q, ©, $5$ and ([Ml) it follows 
that 



A(t)dt = 
o V ' 



-r) 




(60) 



|*> = —(1^(0)0) + 1^2(0)1)) . (65) 



The integral f Q Hi(t)dt in ([55)) can be evaluated by 
performing the integration along the loop shown in Fig. [I] 
with constant velocities. One obtains 



/ H T (t)dt = I 
Jo Jo 



E 2 dt = xTj{0m, 4>m) 



(61) 



with 

i{6m,4>m) 



1 sin 26 M sin 2 4> M + (3 + cos 2(j> M ) 6 M 
3 + 126» 



+ 



sin 2<j>M sin 2 6 M + (3 + cos 29 M ) 4>m 



120 



M 



.(62) 



Here, 6m and 4>m satisfy 4>m(1 — cos9m) = ?7, with r\ the 
solid angle spanned on the sphere during the evolution 
given in (fT5j) . 

From ((601) and (HJ), one may evaluate the RHS of (J58j 
and obtain, in the logical basis {|0), |1)}, the matrix 



-2e l2L ^ L r]S\n ^ 



2e lX 2 7 jysin ^ 



e 2 K , 



, (63) 



Note that the time evolution of |Z?i(0)0) is driven by 
the unitary operator Uo{T) in (TT4"]) . since the interaction 
does not affect the state |0) of qubit B. On the other 
hand, |Z?2(0)1) evolves with the perturbed unitary oper- 
ator U X (T) in (l6"3"l) . Thus, to study the evolution of ^ it 
is useful to introduce the operator 



U(T) = U (T) <g> 



1 




U X (T)® 




1 



(66) 



Here, as usual, Uo(T) and U X (T) act on the qubit A, 
while the other two matrices act on the qubit B and are 
expressed in the |0), |1) basis for B. Then 



W(T)|Di(0)0) = cosr^OO) - sin^lO) 
U(T)\D 2 (0)1) = ^-^(2r/sin^|01)+. 



Let us now evaluate (|26|) at t = T, 



where 

K± = (jl cos ^ ± ixT~f sin ^ 



A* = V V + (xTl) 2 . 

(64) 



p' A = Tr B {P ® PU(T)\^)(y\U*{T)P (g> P) , (68) 



and thus, from |66|) . 



e X 2 7 (2?y cos ?y sin ^ — if_ sin 77) 



(69) 



(2?7 cos r\ sin ^ — i£" + sin ?/) 



The calculation of £ r is now straightforward by using 
(f37|) and ([69]) . Here, we quote the lowest order expansions 
in the coupling x 



( 7X r) 2 sin 4 ( ?? ) 
8t? 2 



+ 0(( X T) 3 )). 



(70) 



Note that the presence of the qubits coupling affects the 
entanglement inducing a quadratically decreasing behav- 
ior in x- Only in the limit 77 — s- (no holonomic trans- 
formation) the entanglement is still preserved, £ r = 1, 
irrespectively on the coupling. Indeed, in this case the 



time evolved state, starting from a maximally entangled 
state (l65l) differs from the initial state by a phase fac- 
tor |*(T)) = (|00)+exp(-ixT)|ll))/\/2, which does not 
affect the entanglement. 

The fidelity (|4"Tj) can now be evaluated for U (T) given 
by ([66]) (see also ([67)) ). One finds 

T = 1 - (7 + 8t7 2 + 8?7sm(2?7) - 8cos(2?y) + cos(4r?)) x 



0((XT) J )) 



(71) 



256?7 2 

Note that, differently from the entanglement, even in the 
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t7 



1.0000 r 



0.9995 



0.9990 



0.9985 



0.9980 [ 



10 



50 



100 



150 



200r 




0.9980 



200r 



FIG. 3: Left: Entanglement as a function of T/t for a = 0.1 and \t — 10 -3 (grey dashed line); \ T ~ 5 10 -4 (black dashed 
line). Right: Fidelity as a function of T/t for a = and \ T ~ 10~ 3 (grey solid line), a = and x r = 5 10~ 4 (black solid line), 
a = 0.1 and x r = 10~ 3 (grey dashed line), a = 0.1 and xi" = 5 10 -4 (black dashed line). In all the plots the loop is as in Fig. Q] 
and with the choosen 9m and <f>M we have 7 = 0.75. 



absence of holonomic transformation (77 = 0) the fidelity 
is affected by the coupling, i.e., 



and the evolution operator which extends 



is 



T\ 



r)=0 



1 



{ixTf 



0((XT) 3 )). 



We shall now take into account the effect of the para- 
metric noise and consider the case of a perturbation to 
the logical NOT operation with rj = tt/2. The perturbed 
unitary operator U x is obtained from (|63[) by means of 
the substitutions 



U{T) = U(T) $ 



1 




U X {T)® 




1 



(72) 



V -> V = 2 + S V 

M -> H = V(7r + 2^) 2 + ( X T 7 ) 2 

T ( ft M 
K± —y K± = {ft cos — ± txT^ sin — 



with i7(T) given in Eq. ([2"2"]). whence 



Pa = 



2/1 



cos (5^ + 2?7 sin (5r/ sin ^ 



A'_ cos (5?y + 2?7 sin <5r/ sin ^ 



/' 



(73) 



Thus, from (f37|). ((4T|) . and l[73]). we finally obtain (at 
the lowest orders in the errors) 



T 



1 - 



1 



(tt - AS V ) { lx Tf 



2tt 3 



(Jr/) 2 



{ixTf 



32ir 2 



(74) 
(75) 



It should be observed that £ r is jointly reduced by cou- 
pling x an d parametric noise Sij (yet £ r = 1 for x = 
and Sr) 7^ 0, in agreement with the results of Sec. HID), and 
that the parametric noise may increase the entanglement 
with respect to the bare case with 77 = 7r/2 (in accordance 
with the behavior discussed in Eq. ([70]) . where £'" — > 1 for 
r/ — > 0). Moreover, note that the fidelity T is influenced 



independently by the two sources of noise. Indeed, in ad- 
dition to the decrease induced by the coupling, already 
accounted by (|7i"|) . the fidelity is also depressed by the 
geometric perturbation Sr/ in agreement with the results 
in Sec. [TTTJ 

The model of parametric noise introduced in Sec. IIVI 
allows for a more quantitative analysis of (|74|) and (|75[) . 
According to this model, 5rj = a^J t/T (cf. Eq l|49p). with 
a 2 the variance and r the time scale of the noise fluctu- 
ations; the condition of small perturbations, Srj -C tt/2, 
implies ay/ t/T <C 7r/2. The behaviors of entanglement 
and fidelity are represented in Fig. [3l more precisely: 

• Fig. [3] (left) shows the entanglement (|74|) as a func- 
tion of the final time T for fixed a and different 
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values of x T - The value of 7, denned in (|62|) . de- 
pends on the loop chosen with 0.66 < 7 < 1. Here, 
7 = 0.75. The behavior suggests a way to pre- 
serve entanglement: choosing not too long evolu- 
tion times the dominant source of error can be min- 
imized. In particular an evolution with T < 200r 
causes an error on the entanglement smaller than 
1.5 10~ 3 . 

• Fig- El (right) shows the fidelity. When a = 0, i.e. 
no parametric error, the fidelity shows a quadrati- 
cally decreasing behavior MS M function of T/r de- 
pending on the values of x T (solid lines). The pres- 
ence of a parametric noise changes qualitatively the 
behavior for small T/r. In fact, in this region in- 
dependently of the x coupling, the fidelity drops 
because the geometric error dominates. By increas- 
ing T/t there is an intermediate region where the 
fidelity increases before fast decreasing, when the 
X error prevails. Intermediate times evolution are 
then the more efficient to preserve the fidelity. A 
good choice is 50r < T < 100t. In this range 
of adiabatic times, both entanglement and fidelity 
have errors of order 5 10~ 4 . 

To conclude, the constraints on the time scales of our 
model that allow to construct holonomic gates which pre- 
serve both entanglement and fidelity and arc consistent 
with the adiabatic approximation are: 

IT 1 < 50r < T < 100r < x" 1 (76) 

with the additional request that a <C \/ r 50n/2. 

VI. CONCLUSIONS 

In the present paper we studied a noisy two-qubit sys- 
tem in order to understand if and how the holonomic 



operators preserve the entanglement. We considered a 
model in which only one of the two qubits undergoes a 
holonomic transformation. Being the two dimensional 
logical space embedded in an extended four dimensional 
Hilbert space we introduced, as possible estimator, the 
reduced logical entanglement which corrcponds to the 
fraction of entanglement in the logical subspace. We also 
calculated the fidelity and compared it with the reduced 
logical entanglement. 

We considered two types of noise: a parametric noise 
that goes off whenever the driving fields are turned 
off, and a coupling noise due to undesired interactions 
between the two qubits. We have shown that the 
holonomic operators are robust under parametric error. 
In particular, the entanglement is preserved under an 
holonomic transformations while the fidelity is affected 
by such a noise but, due to geometric cancellation 
effects, it can reach good values for long times. In the 
presence also of a coupling error, we showed that the 
entanglement is mainly influenced by this noise and 
weakly by the geometric perturbation. Instead, the 
fidelity shows a different dependence on the adiabatic 
time: it is dominated by the parametric noise for not 
too large times and depends on the coupling error at 
larger times. We demonstrated that the intermediate 
time evolutions are the best choice. Within a realistic 
range of physical parameters, we found that for both en- 
tanglement and fidelity the error can be strongly reduced. 
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